Impact of vector-current interactions on the QCD phase diagram 
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Using a nonlocal version of the Polyakov-loop-extended Nambu-Jona-Lasinio model, we investi- 
gate effects of a nonderivative vector-current interaction (relating to the quark-number density) at 
both real and imaginary chemical potentials. This repulsive vector interaction between quarks has 
the following impact on the chiral first-order phase transition: at imaginary chemical potential it 
sharpens the transition at the Roberge- Weiss (RW) end point and moves this critical point toward 
lower temperatures; at real chemical potential, the critical end point moves on a trajectory towards 
larger chemical potentials and lower temperatures with increasing vector coupling strength. The 
conditions are discussed at which the first-order phase transition disappears and turns into a smooth 
crossover. 
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I. INTRODUCTION 

Exploring the phase diagram of quantum chromody- 
namics (QCD) at finite temperature and real chemical 
potential is one of the most interesting and important 
subjects in particle and nuclear physics. Lattice-QCD 
(LQCD) simulations are a powerful method to investigate 
the QCD thermodynamics at zero chemical potential. At 
finite chemical potential, however, LQCD suffers from the 
so-called sign problem which restricts the applicability of 
LQCD to the region of small real chemical potential (/xr) 
and high temperatures (see Ref. [1]). Therefore, model 
calculations (admittedly with substantial ambiguities [2]) 
are used to investigate the phase structure at moderate 
and large /xr. 

A promising strategy for studying the QCD phase di- 
agram at finite /ir is the imaginary-chemical-potential 
matching approach [3]. ft is similar to the usual imag- 
inary chemical potential approach for LQCD [4-7]: 
LQCD data at finite imaginary chemical potential (/Xi) 
are extrapolated to the /xr region by using an analytic 
function. In the imaginary-chemical-potential matching 
approach we extract some important restrictions for the 
model design from the /xi region. This allows us to extend 
the model to the real-chemical-potential region more re- 
alistically (cf. Ref. [8]). The important point is that the 
fj,i region encodes almost all information of the /iR region. 
This fact can be understood through a Fourier transfor- 
mation of the grand-canonical partition function, Z, in 
terms of 9 = [ii/T in the case that the baryon number is 
a good quantum number: 

ZcW = £- f dde- iN * s Z(6). (I) 

Here, Zc(N q ) is the canonical partition function with 
real quark numbers (N q ). Moreover, QCD possesses the 
so-called Roberge- Weiss (RW) periodicity [9] : thermody- 
namical quantities have a periodicity of 2ir/3 along the 
0-axis. This periodicity is described by invariance under 



the extended Z3 symmetry [10] 

9 -> 9 + 2nk/3, $ ->■ e - 27rfe/3 $ , <S^e 27rfc/3 <l>, (2) 

with integer k, where $ is the Polyakov loop and $ its 
conjugate. The RW periodicity enables us to determine 
which interactions are relevant and how strong the cou- 
plings are by comparing model results with LQCD data 
at finite Note, that [i\ can be absorbed in the bound- 
ary angle of the temporal direction of the quark field. 
From this viewpoint, quarks are fermions at 9 — 2nk/N c 
and these become boson-like at 9 = 7r(2fc + l)/A c . There- 
fore, the dual quark condensate was proposed [11] as 
an order parameter for the chiral and the deconfinement 
phase transition. 

At 9 — 7r/3 another characteristic property of QCD, 
the so-called RW transition, arises: this RW transition 
can be related to charge-conjugation (C) or CZ 3 symme- 
try breaking [12]. CZ3 symmetry is explicitly broken at 
finite /xr, but it is not explicitly broken at 9 = 7r/3 be- 
cause of RW periodicity (see, e. g., Ref. [13]). On the RW 
transition line, 0-odd quantities can have a finite value, 
but they vanish for temperatures below the RW end point 
[9]. Therefore, we can interpret 0-odd quantities as order 
parameters of spontaneous C-symmetry breaking. 

The Polyakov-loop-extended Nambu-Jona-Lasinio 
(PNJL) model is a promising approach as it preserves 
the RW periodicity in the same way as QCD. In the 
present study we extend the nonlocal version of the two- 
flavor PNJL model from Refs. [14, 15] by introducing a 
nonderivative vector-current interaction between quarks 
both at imaginary and real chemical potentials. 

This paper is organized as follows: In section II we in- 
troduce the nonlocal PNJL model that is used in our cal- 
culations. In particular, we describe in detail the treat- 
ment of the vector-type interaction in the nonlocal frame- 
work. We show how this approach can be extended to 
imaginary chemical potentials. Section III presents the 
results of our calculations. Section IV closes this work 
with a discussion and a summary. 
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II. FRAMEWORK AND FORMALISM 



Lagrangian density and nonlocality distribution 
functions 



The generic Euclidean action of the two-flavor PNJL 
model is 

S = J dr J d 3 x [q{x){-ifl + mo)q(x) - Ant] 

+ j3VU($[A],$[A};T), (3) 

where q(x) is the two-flavor quark field, Too denotes 
the current quark mass, and D u = d v + iA v = d v + 
ib\ A 4 ' a X a /2 is the color gauge-covariant derivative with 
SU(3) C Gell-Mann matrices A a . The gauge coupling g is 
understood to be absorbed in the definition of A 4 - a . 

The last term in Eq. (3) is the Polyakov-loop-effective 
potential IA 1 multiplied by volume V and inverse temper- 
ature (3 = T -1 , and to be specified later. $ and $ are 
the Polyakov loop and its conjugate, respectively. 

The nonlocal generalization of the PNJL model is char- 
acterized by an interaction featuring nonlocal quark cur- 
rents and densities, as follows [14-18]: 



Ant fa) = G j a (x)j a (x) + J(x)J(x) 

Ja (x) = / d 4 z C(z) q(x + z/2) T a q(x - z/2) 



tf(x) = / d 4 z C{z) q(x + z/2) f q{x - z/2) 



(4) 
(5) 

(6) 



i0 



J(x) = / <Tz F{z) q(x + z/2) — q(x - z/2) . (7) 
J 2k 

The chiral (scalar and pseudoscalar) densities j a (x) with 
a = 0,1,2,3 involve the operators T a = (l,«75-r). The 
overall coupling strength G of dimension [length] 2 is cho- 
sen sufficiently large so that spontaneous chiral symmetry 
breaking and pions as Goldstone bosons emerge properly. 
The term involving the nonlocal quark vector currents 
has a coupling strength G v , again of dimension [length] 2 . 
This G v is treated as a parameter in the present work. 
For orientation, the Fierz transformation of a color-octet 
current-current interaction (induced, e.g., by gluon ex- 
change) gives G/G v = 1/2 (see, e.g., Ref. [8]). 

The term involving J(x) is an additional vector-type 

derivative coupling with q(x') q(x) := q{x'){d ll q){x) — 
(d IM q)(x') q(x) together with a scale k so that the effective 
strength of this term in £ int is G/n 2 . In the following, 
we refer to the interaction induced by J(x) simply as a 
derivative coupling in order to avoid confusion with the 
nonderivative vector interaction which is called vector- 
current interaction from here on. 

The currents Eqs. (5)-(7) include nonlocality distri- 
butions C(z) and F{z). These distributions govern 
the momentum dependences of the quark mass function 



and of the renormalization factor that appears in the 
quark quasi-particle propagator, Z(p 2 )(7 ■ p — M (p 2 )) -1 , 
[8, 15, 17]. The Fourier transform C(p 2 ) of C(z) is related 
to the quasi-particle mass function M(p 2 ) determined by 
the self-consistent gap equation, 



M(p 2 ) = Z(p 2 ) m + aC(p 2 ) 



(8) 



where a is the scalar mean field basically representing 
the chiral condensate (qq). The Fourier transform F(jp 2 ) 
of J^(z) is, in turn, related to the Z factor of quark wave- 
function renormalization, 



Z(p 2 



1 



Hp 2 ) 



(9) 



where v is the mean field induced by J(x) [19]. 

The following four-dimensional momentum-space 
forms of the distribution functions are used in this study: 



C(p 2 ) = J d 4 z exp(-ip-z)C(z) 

( e -p 2 d 2 c /2 (p 2 < A 2) 

" lAf^ (p 2 >\ 2 ) , 
J-(p 2 ) = J d 4 z exp(— ip ■ z)J r (z) 
= exp(-p 2 4/2) . 



(10) 



(11) 



The running QCD coupling a s (p 2 ) determines the 
asymptotic form of C(p 2 ) while its infrared behavior is 
given by a Gaussian parameterization with a character- 
istic length scale dc- The matching of these high- and 
low-momentum representations at an intermediate scale 
A determines the constant J\f. Parameters in both distri- 
bution functions are fitted to LQCD data as described in 
Ref. [151. 



B. Thermodynamical potential 

Consider now the (grand-canonical) thermodynamical 
potential, 



where 



T 



Z = VqVqe 



(12) 



(13) 



is the grand-canonical partition function determined by 
the path integral over the action (3). In the mean-field 
approximation the fields are replaced by their (thermal) 
expectation values. After bosonization, the mean-field 
thermodynamical potential, I^mf °f the nonlocal PNJL 
model, including quark wave-function-renormalization 
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corrections, but in the absence of the the vector-current 
interaction reads 



where 



i=±,0 n= — oo 



r> MF = fii+W($,$;T), 
d 3 p 



(14) 



(2tt)= 



In 



AG 



(15) 



Here cr and v are the mean fields associated with the 
scalar density jo and the derivative vector current J, 
respectively. The first term on the right-hand side of 
Eq. (15) involves the quark quasi-particle energies 



p- 



M 2 (p 2 



(16) 



with dynamically generated masses, Mi(p 2 ) = M(p 2 = 
w ni +P 2 )j determined self-consistently at each shifted 
Matsubara frequency ui rh i with i G {±,0}: 

A 3 A 8 



in ± 



2 v / 3 



in ■ 



V3 



(17) 



are the gauge fields forming the Polyakov loop given 
in Eq. (22). Likewise, the Z factors are understood as 



Z^p 2 ) = Z(p 2 



p ). More explicitly: 



M t (p 2 ) = Z t (p 2 ) m + aC{p 2 = u 2 n ,+p 2 ) 



Zi(p 2 )= \--T(p 2 = 



P 



(18) 
(19) 



At finite temperature T, the Lorentz invariance is bro- 
ken by the thermal medium and the inverse quark quasi- 
particle propagator becomes S~ 1 (p) = —Ai{p)^ip 4 — 
A{p)^iip l + B(p) with A4 7^ A. Here we assume for 
simplicity that the difference between A4 and A is suf- 
ficiently small so that it can be neglected, given that 
the overall influence of wave-function renormalization on 
thermodynamical quantities is not very significant. 

The introduction of the vector-current interaction 
leads to the following modifications of the thermodynam- 
ical potential (14): first, from the bosonization of Eq. (6) 
a quadratic term involving the vector mean field u>, 



UJ 

A~G~ V 



(20) 



is added to the thermodynamical potential (15); sec- 
ond, the chemical potential is shifted according to /1 — > 
fi — C(p 2 = uj 2 ti + p 2 ) uj. The vector mean field ba- 
sically represents the baryon density, g — in the 
form uj = • One important remark is in order: the in- 
dependence does not appear in the distribution functions 
because these functions are introduced in the Lagrangian 
density before taking the mean-field approximation. 



C. Polyakov-loop potential 

In this study, we consider two types of the Polyakov- 
loop effective potentials. The first one is given in 
Ref. [201: 



+ bi{T) ln[l - 6 <£ $ + 4($ 3 + $ 3 ) - 3($ 4>) 2 



(21) 



where $ and $ are represented as 
A\ + AX 



exp 



2T 



exp 



A 3 - A 8 
1 2T 



+ exp 



V3T 



The other one is proposed in Ref. [21] 
W(#,$;T) 



(22) 
(23) 



T 4 



= -bT 



54e r $ $ 



64 (T) ln{l - 6 $ $ + 4($ 3 + $ 3 ) - 3($ $) 2 } 



(24) 



This latter form is obtained from the knowledge of the 
strong-coupling limit of QCD. Recently the details have 
been investigated in Ref. [22, 23]. The nonlocal PNJL 
model with potential (21) is henceforth denoted as model 
A and that with potential (24) as model B. 

It is convenient to introduce a modified Polyakov-loop 
and its conjugate as 



(25) 



as these are RW-periodic quantities. The real and imag- 
inary parts of $ serve as order parameters of the decon- 
finement transition and spontaneous C-symmetry break- 
ing because Im ^> is a 0-odd quantity, just like the quark 
number density. We use Im "J/ as the order parameter of 
C-symmetry breaking. As mentioned in the introduction, 
this 0-odd quantity serves as an exact order parameter 
at 6 = 7r/3 because there C-symmetry is not explicitly 
broken. 



D. Parameter setting 

In the PNJL model the pion mass and its decay con- 
stant are used to fix parameters in the NJL sector of the 
Lagrangian. These parameters are taken from Ref. [15]. 
The vector-current interaction, at mean-field level, has 
no influence on the thermodynamics at fj, = 0. An esti- 
mate of the coupling constant G v can therefore only be 
provided by comparison with (restricted) LQCD infor- 
mation at nonzero chemical potential. For guidance, we 
can use the LQCD value for the ratio Trw/Tc ~ 1.05 [4], 
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where Trw is the critical temperature of the Roberge- 
Weiss end point (at — 7r/3) and T c is the crossover 
temperature at 6 = 0. 

The coefficient functions bi (T) and bi(T) in (24) are 
parameterized such as to reproduce pure-gauge-LQCD 
results (Refs. [15, 16]). When we use T = 270 MeV for 
the confincmcnt-deconfinemcnt transition temperature in 
the pure-gauge case, the resulting crossover transition 
temperature when including quarks in the PNJL model is 
slightly higher than the LQCD prediction. Alternatively, 
we also use T = 240 MeV in order to reproduce T c ~ 
190 MeV, as suggested in Ref. [24]. 

The parameter a in (24) is fitted to reproduce the crit- 
ical temperature in the pure gauge limit and its value 
is a = 664 MeV. The remaining parameter b is defined 
to reproduce the pseudo-critical temperature with dy- 
namical quarks. To reproduce T c ~ 190 MeV from the 
two-flavor LQCD data, we take b — 0.01. 



III. NUMERICAL RESULTS 

Figure 1 displays the T-dependence of the chiral or- 
der parameter and the real part of \IL Here we show 
the results of model A with T = 240 and 270 MeV and 
that of model B with b = 0.01 and 0.02 at /i = 0, in 
both cases with G v = 0. In all cases the transitions 
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FIG. 1. T-dependence of the chiral order parameter and the 
real part of \P at fj, = with G v = 0. The solid and dotted 
lines denote the result of model A with T = 240 and 270 MeV, 
respectively. The left figure shows the result of model B with 
b = 0.01 and 0.02. 

are crossovers. The transition temperatures for the chi- 
ral and deconfinement crossovers almost coincide. This 
property comes from the entanglement of the chiral and 
deconfinement transitions through the nonlocality distri- 
bution functions. 

Figure 2 shows the T-dependence of Im W at 9 — ir/3. 
From the right figure we see that there is a first-order 
RW end point in the case of b = 0.02 which turns into 
second-order for b = 0.01. The LQCD data [25, 26] sug- 
gest that the order of the RW end point is first-order at 
sufficiently small rriQ. From this perspective, b = 0.01 is 
not a suitable choice, but this situation can be modified 
as shown below. 
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FIG. 2. Left (right) figure: To-dependence (fa-dependence) of 
Im as function of temperature T at 9 = it/3 with G v = 0. 
In the left figure, the solid and dotted lines denote the result 
of model A with T = 240 and 270 MeV, respectively. In the 
right figure, the solid and dotted lines denote the result of 
model B with b — 0.01 and 0.02, respectively. 



In the previous figures we have ignored the vector- 
current interaction. As a consequence, the ratio Trw/T; 
exceeds the LQCD prediction [4, 7]. Choosing G v = 
0.4 G in model A with T = 240 MeV, this ratio be- 
comes T RW /T C = 1.04 [8]. Model B with G v = 0.4G 
and b = 0.01 leads to T RW /T C ~ 1.05. These values 
are close to the LQCD result. Figure 3 shows again 
the T-dependence of the imaginary part of the modified 
Polyakov loop (Im $). One observes that the transition 
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FIG. 3. T-dependence of Im * at = tt/3 for G v = 
and 0.4 G. The solid and dotted lines denote the result with 
G v = 0.4 G and G v = 0, respectively. Left figure: result 
of model A with T = 240 MeV; rigth figure: model B with 
6 = 0.01. 

at Trw becomes first-order when introducing the vector- 
current interaction in model B with b = 0.01. Henceforth, 
we only refer to the results of model A as both models 
lead to almost identical results. 

Finally, we study the G v -dependence of the position of 
these critical point in the T-fi phase diagram. In Fig. 4, 
the circles and triangles represent the positions of the 
critical end points for the nonlocal version of the PNJL 
model, respectively. Results are presented for G v /G = 0, 
0.25, 0.4 and 0.45, respectively. In the local PNJL model, 
the critical end point disappears or shifts to very small 
T when considering a realistic range 0.25 < G v /G < 1 
for the vector coupling strength. In the nonlocal PNJL 
model, the location of the critical point has a less pro- 
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nounced dependence on temperature, at least for small 
G v . This behavior can be traced to the weakening of 
the NJL interaction by the nonlocality distribution. The 
downward trajectory of the critical point becomes very 
steep, however, once G v /G reaches values of 0.4 and be- 
yond (see Fig. 4). Around G v /G ~ 0.5, the canoni- 
cal ratio corresponding to an effective quark-quark in- 
teraction induced by color-octet (gluon-exchange) cur- 
rents, the critical points tends to disappear altogether 
and the first-order phase transition turns into a continu- 
ous crossover. 
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FIG. 4. Phase diagram in the ^ir-T plane. The circles and 
triangles represent the position of the critical end point for the 
nonlocal and the local PNJL models, respectively. The white, 
blue, black blue and black symbols correspond to results with 
G v /G = 0, 0.25, 0.4 and 0.45, respectively. 



IV. SUMMARY 

In this study we have investigated the impact of a 
(nonderivative) vector-current interaction in the nonlocal 
PNJL model at real and imaginary chemical potentials. 
The presence of the vector-current interaction makes the 
transition at the Roberge- Weiss end point more pro- 
nounced. The RW end point becomes first-order, con- 
sistent with recent LQCD simulations. 

The location of the critical point in the phase dia- 
gram for real chemical potentials is highly sensitive to 
the vector coupling strength G v . In the nonlocal PNJL 
model used here, the critical point tends to be elimi- 
nated in favor of a continuous crossover once the ratio 
of vector-to-scalar couplings is increased toward and be- 
yond G v /G = 1/2, the value characteristic of an effec- 
tive gluon-exchange interaction between quarks. Qualita- 
tively similar tendencies are found in recent related work 
[27] and in a 2+1-flavor study using the local PNJL model 
[28] in which the disappearance of the chiral first-order 
transition turning to a crossover is indicated already at 
values G v /G < 1/2. 
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